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Abstract 

Recent experiments aimed at probing the dynamics of excitons have revealed that semiconducting 
films composed of disordered molecular subunits, unlike expectations for their perfectly ordered 
counterparts, can exhibit a time-dependent diffusivity in which the effective early time diffusion 
constant is larger than that of the steady state. This observation has led to speculation about 
what role, if any, microscopic disorder may play in enhancing exciton transport properties. In this 
article, we present the results of a model study aimed at addressing this point. Specifically, we 
present a general model, based upon Forster theory, for incoherent exciton diffusion in a material 
composed of independent molecular subunits with static energetic disorder. Energetic disorder 
leads to heterogeneity in molecule-to-molecule transition rates which we demonstrate has two 
important consequences related to exciton transport. First, the distribution of local site-specific 
diffusivity is broadened in a manner that results in a decrease in average exciton diffusivity relative 
to that in a perfectly ordered film. Second, since excitons prefer to make transitions that are 
downhill in energy, the steady state distribution of exciton energies is biased towards low energy 
molecular subunits, those that exhibit reduced diffusivity relative to a perfectly ordered film. 
These effects combine to reduce the net diffusivity in a manner that is time dependent and grows 
more pronounced as disorder is increased. Notably, however, we demonstrate that the presence of 
energetic disorder can give rise to a population of molecular subunits with exciton transfer rates 
exceeding that of subunits in an energetically uniform material. Such enhancements may play an 
important role in processes that are sensitive to molecular-scale fluctuations in exciton density 
field. 
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I. INTRODUCTION 


Excitons are Coulombically bound electron-hole pairs that mediate the transport of en¬ 
ergy in many molecular scale systems. The dynamics of excitons are fundamental to pro¬ 
cesses such as photosynthetic light harvestin^i®, photocurrent generation in solar cell^^®^, 
and photoluminescence in light emitting diodes (LEDs)PlEI. The presence of molecular 
disorder is generally thought to hinder exciton transport propertie^^®^, for instance by re¬ 
ducing diffusivitjff^^, however, it has been speculated that in certain cases disorder might 
serve to enhance exciton transporlP^^. 

The hrst spectroscopic evidence of quantum coherence in photo-excited light harvest¬ 
ing complexe^ inspired a flurry of studies aimed at understanding exciton transport in 
molecular system^®®^. Theoretical model studies using quantum master equatioi J^^hs i 
path-intergral formulatiorP^ESl have indicated that the highly efficient transport of excita¬ 
tion energy in biological and nanoscale systems can be attributed to an effective interplay 
between the two regimes of exciton transport: the local wave-like regime and the non-local 
hopping regime. When electronic coupling between molecular subunits is strong, excitons 
tend to delocalize and their dynamics involve the collective interplay of nuclear and elec¬ 
tronic degrees of freedom. This type of exciton dynamics is referred to as coherent, which 
differs qualitatively from the incoherent dynamics that arise when intermolecular coupling 
is weak in which localized excitons undergo discrete hops between different molecular sub- 
unit^HmSl There are many examples of systems including some organic semiconductors,^^®^ 
colloidal quantum dots,l^®^ and certain molecular aggregated^^® where exciton transport 
is dominated by incoherent-type dynamics. 

In photosynthetic systems, synergy between local excitonic coherence and thermal fluc¬ 
tuation from the environment results in energy transfer efficiency exceeding purely coherent 
or purely incoherent transfer.l^ The relative manifestation of these two regimes is governed 
by static and dynamic disorder .SUZEH The effects of disorder in coherent exciton dynamics 
has been well studied in the context of general theoretical models. While static disorder has 
been shown to arrest diffusion via the phenomenon known as Anderson localization,*^^*^ it 
has also been shown that this effect can be removed by the presence of dynamic disorder, 
for example, in the form of a system-bath couplingj ^^ * ^^ * Relatively little attention has been 
given to what role, if any, disorder plays in enhancing exciton transport in the incoherent 
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regime. 

In this article, we explore the relationship between static disorder and the dynamics of 
excitons in the context of a simple but general model for exciton diffusion in the limit of 
small intermolecular coupling. We demonstrate that although site energetic disorder leads 
to a reduction in the average steady state diffusivity of excitons, microscopic heterogeneity 
gives rise to regions of the material for which exciton mobility is enhanced relative to that in 
the absence of disorder. We go on to show that the heterogeneity in microscopic transport 
properties is correlated with local excitation energy, and that this heterogeneity is responsible 
for the time-dependent diffusivity that has been observed by transient photoluminescene 
spectroscop}^^^® as well as by more recently developed time-resolved optical microscop}^^^®. 
In the next section we describe our model for disordered incoherent exciton dynamics. 


II. MODEL 


We employ a simple and general model of a semiconductor him composed of a collec¬ 
tion of disordered molecular subunits. Forster theory provides the theoretical basis for our 
treatment of incoherent exciton dynamicd^. In our model, excitons diffuse by performing a 
series of discrete hops between intermolecular subunits. Hopping rates are governed by the 
Forster rate equation. 


; ^ ^ f Ro{£d,£a) 

Kda(«, £d,£a) = - 


( 1 ) 


where d is the hopping distance, £d and Ea are the absorption (excitation) energies of the 
donor and acceptor sites, respectively, r is the experimentally determined exciton lifetime, 
and Rq is the Forster radius. The dependence of exciton hopping rate on the excitation 
energies of molecular subunits is contained within the expression of the Forster radius. 
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where n is the refractive index, c is the speed of light, h is the reduced Planck constant, r] 
is the photoluminescence quantum yield, and k is the transition dipole orientation factor. 
The term in the integral contains the overlap between the normalized emission spectrum, 
/{e^Ej)), of the donor molecule with site energy, £d, and the absorption spectrum, a{E]EA), 
of the acceptor molecule with site energy ^a- The line shapes of a{E]EA) and /{e-jEb) 
are taken to be Gaussian with a mean of Ea and ~ respectively, and each with a 
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standard deviation eqnal to the homogeneons line width cxh. For a given molecnlar snbnnit, 
the difference in energy between the absorption and emission peaks is given by the Stokes 
shift, Ass, which reflects the rapid (on the timescales of intermolecular excitonic transitions) 
electronic and nnclear relaxation that follows the excitation. Simulations are carried out in 
the dilute limit of exciton density. 

We incorporate disorder into our model by varying the site energy (he., Sa and Eb in 
Eqs. [^and|^ of individual molecular subunits. For this study the model system consists of 
a two-dimensional array of 2500 hexagonally closed packed molecular subunits with absorp¬ 
tion peak energies drawn randomly from a Gaussian distribution with mean e and standard 
deviation aih, corresponding to the inhomogeneous linewidth. Each molecular subunit is 
assigned a fixed transition dipole vector (used to generate k, in Eq. oriented randomly on 
the surface of a unit sphere. The results presented here utilize parameters and a system ge¬ 
ometry that were selected to correspond to the thin film of CdSe/ZnCdS core-shell quantum 
dots (the subject of recent experiments^. Nonetheless we expect the qualitative features 
of our model to apply generally across systems exhibiting incoherent exciton transport. The 
model system is described in greater detail in Ref. [3H1 


III. LOCAL EXCITON HOPPING RATE DEPENDS PRIMARILY ON DONOR 
SITE ENERGY 

To understand exciton transport from a microscopic perspective, we consider the phe¬ 
nomenological rate for exciton transfer from an individual electronically excited molecular 
subunit. For condensed phase systems which exhibit incoherent exciton transport, this rate is 
dominated almost entirely by the sum of nearest neighbor molecule-to-molecule transitions, 
which are uniform and constant in a perfectly ordered film but vary upon the introduction 
of disorder. In our model such variations arise from the dependence of Eqs. [T] and on 
energetic disorder, through the variables £d, £a, and spatial disorder, through the variables 
d and k. From the perspective of an electronically excited donor molecule, the local envi¬ 
ronment provides a sampling of the system heterogeneity. For multidimensional systems, 
the dependence of exciton transfer rate on the environmental variables Ea, d, and k are mit¬ 
igated by the population of many nearest neighbors. The qualitative effect of site energetic 
disorder on exciton transport can therefore be understood primarily in terms of the donor 
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site energy by considering the average molecule-to-molecule transfer rate from a molecular 
subunit with site energy to a nearest neighbor site, given by 
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where do is the nearest neighbor distance and P{£a) is the distribution of site energies in 
the system, 
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which is in this case Gaussian with mean e and variance (Jih- It turns out that even in 
two dimensions the physics contained within Eq. are sufficient to understand the general 
qualitative effects of energetic disorder on exciton transport. 

For systems with symmetric and singly peaked energetic polydispersity, as is the case 
in this study, the expression in Eq. can be simplified in the context of a mean field 
approximation to take the form fcMF(£D) = fcDA(do,£ d, k). Figure contains a plot of 
^mf(£d), where it is shown that the mean exciton hopping rate is a non-monotonic function 
of donor site energy. The maximum transition rate is achieved when the donor site energy 
is above the average, specifically when = ^ + where the magnitude of the donor- 

acceptor spectral overlap integral in Eq. |^is optimized. (See Fig. &•) 

The average hopping rate of a material reflects the distribution of local hopping rates. 
In the next section we evaluate the system-wide average hopping rate and discuss how it 
depends on the details of static disorder. 


IV. THE EFFECT OF DISORDER ON THE MATERIAL-WIDE AVERAGE EX¬ 
CITON HOPPING RATE 

The net effect of site energetic disorder on the material-wide exciton hopping rates can be 
evaluated by comparing films with and without energetic disorder. From Eqs. and we 
see that in the absence of energetic disorder (he., P{£a) = ^{ea — £)), the nearest-neighbor 
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FIG. 1. (a) The average pairwise energy transfer rates, ^mfCed) between a donor molecule with 
site energy ed and an acceptor molecule with average site energy e. (b) The spectral overlap of 
the donor (red) and the acceptor (blue) subunits for selected values, numbered (I) through (IV), of 
^mf(£d) as shown in (a). The difference between the acceptor site energy (same as its absorption 
energy) and the donor emission energy (site energy offset by the Stokes shift, i.e., ed — Ass) is given 
by As = e — sd + Ass- Subunits with site energies between e and e + 2Ass, indicated by the shaded 
green region, exhibit exciton transfer rates that are enhanced relative to a perfectly uniform film. 
Hopping rates are calculated using Eq. with cjih = O.StTh, Uh = 42 meV, Agg = 38 meV, n = 1.7, 
r] = 1, do = 8 nm, = (k^) = 2/3, and r = 10 ns. 
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exciton transfer rate is constant and is equal to 
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which includes the absorption cross section a at d, cT(e), and the average dipole orientation 
factor, {k^) = 2/3. The expression in Eq. [^reflects the the assumption that the term 
does not vary significantly over the Gaussian term with respect to e in Eq. jg 

valid in cases where the linewidth is smaller than the band gap. 

For a him with energetic disorder, the nearest-neighbor exciton transfer rate depends on 
exciton energy, as illustrated in Fig. [pi. Notably, the molecular subunits with site energies 
between e and £-l-2Agg exhibit exciton transfer rates that are enhanced relative to a perfectly 
ordered him. However, a consequence of the non-monotonic turnover in exciton transfer 
rates for > ^ + Agg is that the system-wide average exciton transfer rate is expected to 
be reduced relative to that of an energetically ordered system. More formally, the hlm-wide 
average exciton transfer rate is given by 
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Since and P^Sa) are Gaussian density of states centered at e for both donors and 

acceptors as dehned in Eq. jp Eq. [^ can be rewritten as 


{k) = 


X exp 



(9) 


To evaluate Eq. [p we hrst integrate with respect to Sa and Sb and then e: because the terms 
involving ea and £d are independent from one another while e is dependent on other two 
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variables. This leads to 
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where {3 = cXih/cTh is the ratio of static disorder to the dynamic disorder. Again assnming 
that the term does not vary signihcantly over the Ganssian term with respect to e in 


Eq. we can simplify the above eqnation as 
(fc) = 




which can be rewritten in the terms of the original variables as 
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By comparing Eqs. and [T^ we hnd that the average exciton transfer rate is decreased by 
disorder (he., {k) < /cDA(do,A^)) if 
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The right hand side of Eq. 13 is a monotonically increasing fnnction of /d^, and also is greater 
than or eqnal to 2, with the eqnality occnrring when /3 = 0. Therefore, for any non-zero 
valne of f3, the average exciton hopping rate for an energetically disordered him can only be 
enhanced relative to that of an ordered him {e.g., (k) > /cda {dQ,e,e)) when 

—V>2. (14) 

CTh / 


The valne on the right hand side of Eq. 13 grows with f3 and is thns larger than 2 for typical 


systems {e.g., 2.5 for the system considered here in which (3 = 0.8). Eqnations 13 and 14 tell 
ns that energetic disorder only serves to enhance the hlm-wide average exciton transfer rate 
when the ordered system has poor donor-acceptor spectral overlap, in which case disorder 
serves to facilitate energetic resonance between some nearest neighbor pairs. An analogons 
type of analysis can be carried ont with respect to spatial disorder in which the molecule- 
to-molecnle transition rates are modnlated throngh the d-dependence (or ^-dependence) of 
Eq. The resnlts are analogous to those described above; the presence of spatial disorder 
typically leads to a reduction in the hlm-wide average exciton transfer rate. 

The observed rate of exciton transfer is only equal to the material-wide average when 
excitons are are uniformly distributed across the ensemble of molecular subunits. Energetic 


























disorder, however, has another important effect on exciton transport in that it facilitates the 
time-varying dissipation of excitation energy, thereby biasing excitons towards lower-energy 
molecular subunits, those with reduced effective hopping rates. (See Fig. In the next 
section we employ Monte Carlo simulations to explore time-dependent effect of energetic 
relaxation on observed exciton diffusivity. 


V. TRANSIENT ENERGETIC RELAXATION LEADS TO TIME DEPENDENT 
EXCITON DIFFUSIVITY 

Each energy transfer event is associated with small energetic losses. In our model, the 
dissipative mechanism appears in the form of the Stokes shift, Agg, designed to mimic the 
rapid internal electronic and nuclear relaxation that follows molecular excitation. Because 
of the Stokes shift, the most probable transitions—those which maximize overlap between 
donor emission and acceptor absorption line shapes—are to acceptor molecules with exci¬ 
tation energies that are Agg lower than that of the donor molecule. Therefore, the initial 
diffusion of randomly populated excitons is accompanied by a decrease in average exciton 
energy. 

Using a standard kinetic Monte Carlo algorithm we have studied the dynamics of excitons 
by analyzing the trajectories of individual non-interacting excitons, with a particular focus 
on the relationship between exciton energy and diffusivity. Figure compares the time 
dependent average displacement of a uniform distribution of excitons in an energetically 
disordered film to that in a film without energetic disorder. We observe that excitons within 
the disordered film (see the curve labeled as ‘disorder’ in Fig. [^) exhibit reduced diffusivity 
when compared to a perfectly ordered film (see the curve labeled ‘uniform’ in Fig. |^). In 
addition, we show the average displacement for the subpopulation of excitons initialized on 
high energy molecular subunits, those with excitation energies Si = e + Agg, or low energy 
molecular subunits, those with excitation energies = e — O.SAgg. We observe that the 
initial behavior of either subpopulation is consistent with the results plotted in Fig. 
namely that excitons with low initial energies diffuse more slowly than the ensemble average 
while excitons with high initial energies diffuse more rapidly than the ensemble average. 
The high energy excitons even manifest a short time enhancement in diffusivity relative to 
that of an energetically uniform film. The enhancement in diffusion for high-energy excitons 
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is short lived because this population of excitons rapidly thermalizes, preferring to occupy 
lower energy molecular subunits whose effective hopping rates no longer reflect their initial 
enhancement. (See Fig. & .) Consequently, over timescales that permit multiple site-to- 
site transitions, the average distance traveled by excitons, regardless of their initial energy, 
becomes less than the average distance traveled by excitons in the absence of energetic 
disorder. 

Evidence for the rapid dissipation of energy for the diffusion-enhanced excitons can be 
seen by analyzing the transient energetics of simulated exciton trajectories (represented by 
line color in Fig. |^). Specihcally, the line color in Fig. corresponds to the mean site 
energy, averaged over a given population of excitons, at a specihc time t. In contrast to the 
rapid energetic decrease observed in the population of excitons with initially high energies, 
excitons with initially low energies exhibit an initial increase in energy, indicating that the 
lowest energy sites are not traps but are dynamically repopulated. The transient energetics 
can be understood more deeply by considering the time-dependent distribution of exciton 
energies (see Fig. &)■ 

In Figure]^ we plot P{e), the probability that an exciton resides on a molecular subunit 
with excitation energy e, resolved as a function of time, where at t = 0 excitons are uniformly 
distributed across the all molecular subunits. We observe a pronounced transient red-shift 
starting from an initially uniform distribution (black line in Fig. [^) to a hnal near steady 
state distribution (red lines in Fig. [^), the result of individual excitons performing energy 
lowering molecule-to-molecule transitions. During the non-equilibrium relaxation, the pop¬ 
ulation of excitons residing on molecular subunits with enhanced hopping rates (relative 
to systems without energetic disorder) decrease rapidly with time, which is accompanied 
by a corresponding increase in the population of excitons that reside on sites with reduced 
hopping rates (relative to a system without energetic disorder). The combined effect is a 
time-dependent lowering of net exciton diffusivity as more excitons occupy molecular sub¬ 
units with low energy and reduced inter-molecular transition rates. 

The relationship between static energetic disorder and transient exciton energetics has 
been previously explored in detail by Basslei®. The phenomenon is experimentally mea¬ 
surable through the time-dependence of the fluorescence spectrumP^ESl Similarly, the rela¬ 
tionship between static energy disorder and exciton transport is experimentally accessible 
through the time- and position-dependence of the fluorescence spectrumpSEoEl -v^here the 
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FIG. 2. Exciton dynamics within energetically disordered molecular film, (a) The average diffusion 
length of excitons in a disordered film plotted relative to those in an energetically ordered him, 
the latter denote by (|i?|)Q(t). The four curves correspond to the ordered him result, an average 
over excitons initialized uniformly across within a disordered him (the curve labeled disorder) and 
an average over the subpopulation of excitons initialized on either high energy (labeled high e*) 
or low energy (labeled low Cj) sites. This quantity has been normalized by the diffusion length in 
the ordered him measured at the exciton’s lifetime ((|i?|)Q(t = r)). The average site energy visited 
by the exciton is indicated by its correpsonding color. (Blue is high and red is low energy.) (b) 
Time-evolution of the probability density distribution of excitons with energy e. Time points are 
indicated by its corresponding color. (Black is the initial state, and red is the state in which most 
excitons have already decayed.) 
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FIG. 3. Mean square displacements (MSD), of excitons in molecular films with various degrees 
of site energetic disorder. MSD is normalized by the square of the nearest neighbor distance, do, 
such that {B?)/do^ is proportional to the square of the number of exciton hops. In the presence of 
static disorder, is shorter than that of the uniform energy case, cxih = 0, at all times. 

transient lowering of net exciton diffusivity manifests itself in the mean squared displace¬ 
ment (MSD), In the absence of disorder (he., (Tih = 0), exciton transport is diffusive 

(he., {R^) oc t) at all times.®' On the other hand, static disorder gives rise to subdiffusive 
transporlpSE3B3 (he., {R^) oc where a < 1) with a MSD that grows non-linearly at short 
times and only becomes linear at longer times when the system has reached a dynamic 
equilibrium. As shown in Fig. the net effect of disorder is to reduce the mean squared 
displacement of excitons, especially at long times. Additionally, the extent of initial devia¬ 
tions from the = 0 case depends sensitively on the extent of static disorder, as does the 
steady state diffusion constant (he., the slope of at long times). 


VI. CONCLUSIONS 

In summary, using a simple Forster-type model of incoherent exciton dynamics in an en¬ 
ergetically disordered array of molecular subunits, we have shown that site-energetic disorder 
leads to variations in transport properties of individual molecular subunits. The variations 
correlate with site energy, which can be understood in terms of spectral overlap integrals, 
and have the net effect of reducing the average molecule-to-molecule transition rate relative 
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to that of a perfectly ordered material. Furthermore, in the presence of disorder, the steady 
state distribution of exciton energies is not identical to the distribution of site energies due 
to the dissipation of electronic energy that accompanies typical molecule-to-molecule tran¬ 
sitions. Consequently, the steady state distribution of excitons is biased towards sites with 
lower than average excitation energies, and thus with exciton transfer rates that are reduced 
relative to that in an ordered material. Hence the observed exciton diffusivity is further 
lowered by transient exciton diffusion. The answer to the the question posed in the title 
of this article is length and time scale dependent. For macroscopic systems the answer is 
no: the average effect of disorder is to reduce overall exciton transport. However, at the 
nanoscale the answer is yes: disorder can serve to enhance exciton mobility in certain micro¬ 
scopic regions. Such local enhancement may have played a role in optimizing the structure 
of chromophore network in light harvesting complexes.^ Considering such microscopic en¬ 
hancements (along with their diminishing counterparts) will become increasingly important 
as devices continue to incorporate nanotechnological elements. 
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